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Abstract. - The Continuous Time Random Walk (CTRW) formahsm is used to model the non- 
Poisson relaxation of a system response to perturbation. Two mechanisms to perturb the system 
are analyzed: a first in which the perturbation, seen as a potential gradient, simply introduces a 
bias in the hopping probability of the walker from on site to the other but leaves unchanged the 
occurrence times of the attempted jumps ( "events" ) and a second in which the occurrence times 
of the events are perturbed. The system response is calculated analytically in both cases in a non- 
ergodic condition, i.e. for a diverging first moment in time. Two different Fluctuation-Dissipation 
Theorems (FDTs) , one for each kind of mechanism, are derived and discussed. 
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i Introduction. — Complex physical systems provide a 
ilumber of challenges that cannot be addressed with the 
traditional approaches of equilibrium statistical physics. 
One phenomenon that reveals a striking departure from 
traditional statistical mechanics is Blinking Quantum 
Dots (BQD) [1]. BQD systems jump between two states, 
".light on" and " light off' , thereby being a realization of di- 
dhotomous fluctuations and this two-state model of BQD 
proves to be sufficient to satisfactorily explain the observed 
data. There is consensus that these fluctuations are simul- 
taneously non-Poissonian and renewal [2,3]. This prop- 
erty makes them a paradigmatic example of weak ergodic- 
ity breakdown [4]. The Kubo- Anderson spectral diffusion 
theory [5] is violated [6], and a new theoretical approach 
has to be created to study the emission and absorption 
processes [7]. 

Note that the theoretical approach behind the work of 
Refs. [4,6,7] is inspired to the Continuous Time Random 
Walk (CTRW) theory [8]. However, it is not yet clear if 
alternative approaches based on the adoption of a Hamil- 
tonian prescription are possible. Under the simplifying 
assumption that the "light on" and the "light off' states 
have the same waiting time distribution 'ip{t) the theoret- 



ical CTRW prescriptions are equivalent to the following 
Generahzed Master Equation (GME): 

-p{t)^-l <^{t-t')K.p{t')dt', (1) 



dt' 



where p(t) is a two-dimensional vector, whose components 
p+{t) and P-{t) denote the probability for the system to 
be in the "on" and "off' state, respectively. The matrix 
K is defined by 



K = 



1 -1 
-1 1 



(2) 
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and the memory kernel function <I>(t) is connected to 
ilj{t) by means of their Laplace Transform as: <t{u) = 
u tlj{u)/{l—ip{u)). The formal equivalence between CTRW 
and GME was originally established by Kenkre et al. 
[9] and extended to include renewal aging by AUegrini, 
Aquino et al. [10,11]. 

Is the GME-CTRW equivalence substantial as well as 
merely formal? This is a question of great importance 
for statistical physics. It is well known [12] that starting 
from a Hamiltonian description of a system plus environ- 
ment, that a Liouville equation formalism can be used to 
project the dynamics onto that of the system alone. This 
projected representation is equivalent to the GME. 
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We are convinced that the question of the physical 
equivalence between GME and CTRW is still not yet set- 
tled, and in this letter, by deriving an exact linear re- 
sponse result, using only CTRW arguments, a step to- 
wards solving this controversy is taken. We notice in fact 
that the Fluctuation-Dissipation Theorem (FDT) of the 
first kind [13], which has a Hamiltonian origin, can be 
generalized [14], within an approach compatible with a 
Hamiltonian origin, to a non-stationary condition with the 
following expression: 



RAs(t,t') ^ f3—{Ait)B{t')) t>t', 



dt 



(3) 



connecting the response function RABit^t') to the corre- 
lation function {A{t)B{t')). 

In Eq.(3) A denotes the variable of interest, B is the 
variable of the system through which the external pertur- 
bation affects its dynamics and /3 the inverse temperature. 

To adapt this theoretical prediction to the BQD process, 
we assume A = B = £,{t), with ^ = -1-1 and ^ = — 1 in the 
state |-f) and |— ), respectively and /? = 1. Consequently, 
the FDT attains the form 



where 



R{t,t') = —C{t,t'), 



c{t,t')^{at)m)- 



(4) 



(5) 



It is important to remark that the response of the two- 
state system, compatible with the CTRW approach, has 
been studied by Sokolov [15] as well as Sokolov and Klafter 
[16], with a treatment that, as we shall see, leads to the 
adoption of Eq. (4). Thus, we may conclude that the 
formal equivalence between CTRW and GME indeed cor- 
responds to physical equivalence. However, we notice that 
the theoretical predictions of these authors all rest on the 
same assumption, i.e. that the perturbation seen as a po- 
tential gradient, simply introduces a bias in the hopping 
probability of the walker in moving between the sites while 
leaving unchanged the occurrence times of the attempted 
jumps. We name this scheme the "phenomenological ap- 
proach" . 

On the other hand, the experimental observation of the 
BQD phenomenon affords information on iplt), which, as 
we shall see, depends on two parameters, ^ and T. The 
power-law parameter ^ is the complexity index and de- 
pends on the cooperative nature of the system. The pa- 
rameter T determines the time scale at which the waiting 
time distribution tjj{t) can be considered to be equivalent 
to an inverse power law with index /i. It is realistic to 
expect that the external perturbation has the effect of in- 
fluencing these two parameters, and especially the param- 
eter T, given the fact that /z has a cooperative origin, and 
consequently is expected to be the more robust of the two. 
To calculate the system's response to perturbation in this 
scheme is a difflcult problem, and the aim of the present 
article is to describe a rigorous way to solve it, within the 



spirit of CTRW. In this letter by calculating analytically 
the response to a perturbation of '0(i) we obtain as main 
consequence the new FDT relation 



R{t,t') 



d 
di 



C{t,t'). 



(6) 



Response and FDT within the Phenomenologi- 
cal Approach. — Let us consider the process of a BQD 
switching between "on" [\+)) and "ofF ([-)) states. This 
process can be modeled as a random walk on a two sites 
lattice. At random times determined by the distribution 
■0(i) the walker attempts a jump to the other site, we 
call these attempts "events". Whenever an event occurs 
the walker can either jump to the other site or remain 
in the same site with equal probability 1/2. Within the 
"phenomenological approach" the perturbation f{t) sim- 
ply introduces a bias in these probabilities while leaving 
the statistics of the events, i.e. '0(t), unchanged. 

Thus, under the influence of the perturbation f(t) the 
walker jumps to the other state or remains in the same 
state with probabilities respectively: 



w±^±{t) = Pr{± ^ ±) = 



1 ± ef{t) 

2 
1 T ef{t) 



(7) 



As mentioned above, this case has already been solved in 
a GME approach in Ref. [15] and within a different ap- 
proach in Ref. [17]. To the end of illustrating our method 
and showing its general validity we rederive the response 
in this case with a CTRW approach. Consider a random 
walk on a two-site lattice which starts at time t = with 
equal occupation of the two sites. The occupation vec- 
tor is denoted by p and the distribution of waiting times 
for the walker between the events by V'(0- Wc make the 
assumption that: 

'/^(T) = ^(l + r/T)-^ (8) 

having in mind systems such as the BQDs whose blinking 
behavior is characterized by a residence times distribution 
which is a power law with diverging first moment (i.e. 
1 < ^ < 2). The occupation at time t will be determined 
by all the possible ways the walker can move in the interval 
of time [0, t] between the sites. In the presence of the 
perturbation, turned on at time t = 0, a bias will occur. 
Let us divide all the trajectories in four groups that we 
indicate with gij , determined by the starting site i at time 
i = and the arrival site j at time t, with i,j = ±. The 
fraction of trajectories gij starting in site i at time t = 
and switching to site j at time t is given by Aij(i)p,;(0), 
with: 

A^oii) ^ H H dh^{h)w,,kAti) ■ ■ ■ (9) 

t 

dtnlpitn - t„-l)wfc„ j(i„)*(t - tn). 



Title 



The trajectories are then sequences of "laminar phases" 
between events in which the walker remains in the same 
site. The i/j{tj — tj^i)dtj are therefore the probabilities 
of occurrence of an event at time t = tj without implying 
necessarily a change of site, while ^(t) is just the survival 
probability, i.e. the probability of occurrence of no event 
in an interval of time of length t. 

Let us focus on the group g++: the trajectories belong- 
ing to this group, at time t, are again in the state |+), so 
they contribute only to the population of the first site at 
time t. Eqs. (7) show that, the first-order contribution 
to the response is obtained by keeping in the terms of the 
sum of Eq. (9) only one perturbed weight 'Wki,ki+iiti+i) 
and leaving the unperturbed weight 1/2 in the remain- 
ing factors . This implies that for the trajectories in the 
group g_|_+ characterized by the same set of n events, only 
the terms with the perturbed jump probabilities before 
the last laminar phase may give an overall nonzero con- 
tribution. In fact, consider trajectories with the same set 
of events (in the group g++) , for every trajectory with a 
laminar phase weighted with the perturbed jumping prob- 
ability of Eq.(7), there will be one trajectory identical in 
every respect apart from the following laminar phase being 
of the same length but of opposite sign (i.e. the walker is 
in the other site). These two trajectories according to Eq. 
(7), are weighted with a jumping probability with opposite 
sign and therefore their total contribution to first order in 
e is null. The only surviving contribution to first order 
in £ will be the one with perturbed jumping probability 
before the last laminar phase, the latter being fixed by 
the group considered. The argument extends to a generic 
group gij , the total contribution of this group to the pop- 
ulation of site j at time t is therefore given by Aij{t)pi{0) 
with Aij (t) = A'^j (t) -|- A^j (i) , that is an unperturbed plus 
a perturbed term. For the response we are interested only 
in the contribution of the perturbed part which at the first 
order, as explained above, simplifies to -~jA']^{t)pi{0) with 



Alit) 



V / diiV'(ii)' 
1 Jo 



(10) 



n=l 



t„-i 



dtniitn - i„_i)-/(i„)*(i - i„) 



i.e. the contribution obtained by keeping a perturbed 
jumping probability only before the last laminar phase. 
Summarizing, the contribution of all the four groups of 
trajectories to the population vector at time t is: 



9++ - 


> -Am 


9+- ~ 


> Alit) 


g-+ - 


> -Alit) 


9— - 


> Alit) 



P+(0) 




P+(0) 

P-(0) 




p^iO) 



(11) 



The response of the system is obtained by calculating 
the difference S of population between the two sites. From 
(11) one derives: 



I](t)=p_(t)-p+(f)=2Ai(i). 
Using Eq. (10) it follows: 

oo „t 

S](f)=£^i?e[/ dt„/(t„)V«(tn)*(t 



(12) 



tn)i (13) 



where ipnitn) is just the n-times convolution of tp{t) and 
the sum runs from n = 1 since the n = term, having no 
events, does not contain any term linear in e. From Eq. 
(13) the response function turns out to be: 

oo 

R{t, t') = J2 Mt')'^it - t') = P{t')'i'{t - t'). (14) 

The result of [17] and [15] is then recovered. The autocor- 
relation function for the unperturbed case is known to be 
equal to the function ^(f,t'), see [11], that is: 



c(t,t') = *(i,0 



dx^{x, t') 



(15) 



where ^(f, t') is the probability that, for fixed t' no events 
occurs until time t > t\ and 



ij{t,t') ^ ip{t) + dxij{t - x)P{x) (16) 

Jo 

is the probability distribution that, for fixed t' , the first 
next event occurs at time t. It follows that 



dt 
From Eq. (16) 



-C{t,t') 



dx—ip{x,t'). 



dip{t,t') 
dt' 



dt 



= 4'{t - t')P{t') 



and then inserting this into Eq. (17) one obtains: 
±C{t,t')^R{t,t'), 



(17) 



(18) 



(19) 



that is the general FDT derivable for a system perturbed 
from a canonical equilibrium. 

Response and FDT in the dynamic approach. 

Let us now consider a scheme in which the bias is due 
to a perturbation of the statistics of occurrence times of 
events. We name this approach "dynamic" because we 
have in mind a dynamic model generating the waiting time 
distribution of Eq. (8) through the equation of motion 



y = "0 2/' 



(20) 



with Qfo = (m ^ l)/2^7 describing a particle moving on the 
interval / = (0, 1] that, everytime it gets to 1, is re-injected 
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back inside the interval / with a random initial condition. 
The arrival of the particle at the border generates the 
"events" that lead the random walk on the two-sites lat- 
tice. As explained in Ref. [17] the perturbation in this 
scheme changes the parameter T and therefore ag in the 
following way: 



a±W = «o(l±f/(t)), 



(21) 



where the sign ± depends on the site at which the walker 
is residing. We remark that in this scheme, when an event 
occurs, the walker jumps to the other site or remains in the 
same site with unchanged probability equal to 1/2. The 
perturbation therefore modifies the unperturbed statistics 
of the events, i.e. ip{t) as given by Eq. (8), and introduces 
an "age dependence" in the statistics of the events break- 
ing the time translation invariance of the renewal process. 
At first order in the perturbation (see Ref. [17] ), V'(0 i^ 
changed into the following conditional probability density: 



^±(r|t') 



VXT)[l±e/(<' + T)] 



(22) 



T 



-M 



T 



-V^M+iM 



t +T 



dxf{x), 



i.e. ^p^{T\t')dT is the probability that, after an event has 
occurred at time t', the following event occurs at t = t' + t. 
-0^+1 (t) is the same as Eq. (8) but with index ^ -f 1. We 
consider the case of harmonic perturbation, then Eq. (22) 
turns into: 



^±(r|t') 



T 



^P{t)[1 ± ecos{uj{t' + t))] 
.M-1 



(23) 



Tio 



'ipf^^i{T)[sin{Lu{t' + t)) — sin{Lut')]. 



Keeping only the first term in Eq. (23) is equivalent to an 
unperturbed distribution of the events followed by a bias 
in the choice of the direction of motion similar to the phe- 
nomenological approach. The second term, besides assur- 
ing normalization breaks this equivalence entangling the 
perturbation with the dynamics. By adopting the same 
notation as in the phenomenological case, the contribu- 
tion of the group gij to the population of the site j at 
time i, is Aij(t)pi(0) with 



00 „t 



dii^^(ti|0)i. 

n=0 ki,...k„ 
t 

t„-i ^ 



1 



where: 



/oo 
dxi)f{x-tn\tn) 



(24) 



is the conditional survival probability of remaining in the 
state ±, with no events in the interval of time i — 1„, after 
an event occurred at time t„. Let us rewrite the survival 
probability in the following way 

*±(i - i„|<„) = *(t - t„) ± e^e{t - tn\tn) (25) 



with 

*e(i - tr^K) = ^^{t- i„|i„) - *'(t - t„|t„), (26) 

where, as can be deduced from Eqs. (24) and (23), 
'^{t — tn) is just the unperturbed survival probability while 
'^^{t — tn\tn) and ^^^(t — <„|i„) are the conditional survival 
probabilities originating respectively from the first and the 
second term linear in e in Eq. (23), that is: 



'i>^{t-tn\tn)=Re[ dxij{x-tn)e 



(27) 



and 



*'(t-i„|t„) 



Tuj 



Im[\ dxij}fi,+i{x-tn){e 



lUJX _ iujtr, 



)] 

(28) 

Adopting analogous considerations as those used for the 
phenomenological case, it can be inferred that the first 
order contribution to the population at time t is obtained 
considering only one laminar phase with the perturbed 
distribution of Eq. (23) while for the others keeping the 
unperturbed one. Also in this case one easily realizes that 
the contribution to the response comes only from those 
trajectories with perturbed last laminar phase, that is by 
replacing all the ^'^(^fc+i — ik\tk) in Eq. (24) with the 
unperturbed probability tp(tk+i —tk) and considering only 
the contribution linear in e coming from ^f{t — t'\t'). In 
this case therefore, the trajectories with no events between 
and t will also give a contribution to first order in e. Let 
us consider first only the contribution of the trajectories 
with at least one event, this contribution is the same for 
both the groups gu and gij and amounts to jA'j^{t)pi{0), 
with: 

^lW = 2 E / dtnMtn)^e{t - t„|tn), (29) 



which applies to the probability vector as: 



9+- 



9- 






P+(0) 




P+(0) 

P-iO) 




P-iO) 



(30) 



Then one has to add the contribution of the trajectories 
with no events which gives: 



*e(0|i)p+(0) 

-*e(0|i)p_(0) 



(31) 



In the end calculating p-{t) — p+{t) one obtains: 

I](i) = 2^1(i) + *,(0|t) (32) 



e / dt'[6{t') + Y,Mt')]'^e{t-t'\t') 

•^0 «=1 
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which is just the real part of the Laplace transform of 
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Fig. 1: Numerical simulation (dots) for E(f) in response to 
a constant perturbation as compared to tiie inverse Laplace 
transform of the analytical expression Eq. (36) in the limit 
a; ^ 0. From top to bottom the curves refer to ip{t) ~ ^^f^ (1 + 
t/ry with ^ = 1.45, T = 1.6 and ^ = 1.25, T = 1.6 (e = 0.1). 



Let us consider the contribution coming only from the first 
term in Eq. (26): the Laplace transform of its contribution 
to Eq. (32) is easily calculated: 

t,is)^eRe[U ^<-;-^-^^\-;-A ], (33) 

which for a constant perturbation, i.e. in the limit uj —^ 0, 
gives the asymptotic value of epsilon. The contribution to 
Eq. (32) of the second term in Eq. (26) amounts to 



S2(S) 



-e Im 

Tuj 



1 



3(1 — ip{s — iuj)) 



~ip^,+i{s - iuj) - (1 - '0^+1 (s)) 



V',,+i(-iw) 



(34) 



Considering that ■!/'^_(.i(s) = l+Ts'ip{s)/{l~^i) and making 
a convenient simplification we get: 



£2(5) = e/jTi 



{iiu)'ip{—iLu) + (s — iLo)ij}{s — iuj) — STp{s) 



-eRe 



—suj{l — ip{s — iuj)) 
'ip{—iuj) — ip(s — iuj) 1/1(5) — -0(5 — iuj) 



s(l — %l){s — itu)) iu){l — ij}(s — iuj)) 



.(35) 



Taking the limit for a; — > carefully so as to obtain the 
case of a constant perturbation, and then taking the limit 
s — > (i.e. t — > cxd), in Eq. (35), the final contribution 
amounts to e(— /x). This contribution added to that of 
the first term given by (33) gives the value e(l — /i). We 
confirm this result with a numerical simulation for the case 
of constant perturbation, the good agreement is shown in 
Fig. 1. 
Remarkably, summing (35) and (33) and simplifying 

t(.).I^,(.) + i^.(.).ei^e[ifc-i^^^], (36) 



'0 



"'-"^''dt\ 



(37) 



where 'ip{t^t') is given by Eq. (16). This expression 
uniquely identifies Tp(t,t') as the response function and is 
straightforwardly extended to a generic perturbation f{t), 
simply by using its Fourier representation and following 
the same steps used to prove Eq. (37). 

FDT in the dynamic case. We conclude that in the 
dynamic case the response to an external perturbation f{t) 
can be written in the form 

ft 



E(t) 



dt'R{t,t')f{t') 



(38) 



with R{t,t') = ■4'{t,t'). It can be seen, by differentiating 
Eq. (15) respect to t, that in this case a different Fluctu- 
ation Dissipation Relation (FDR) is obtained, namely: 

^^C{t,t') = -R{t,t'). (39) 

Of course, if a stationary condition is reached, time trans- 
lation invariance is fulfilled and the two FDRs are equiv- 
alent. In order to address the discussion of how general 
the new FDT may be, we point out that the same result 
has been recently obtained by AUegrini at al. [18] using 
a dynamical model different from that of Eqs. (20) and 
(21), fitting only the constraint of yielding the distribution 
of Eq. (8). On the basis of this we argue that Eq. (39) 
might be indeed a universal result, requiring only that the 
perturbation changes slightly the occurrence times of the 
events of the system, thereby producing a small change of 
the parameter T of Eq. (8). 

Summarizing, in this letter we analyzed the problem 
of the response to perturbation of systems with slow re- 
laxation properties described by non-Poissonian renewal 
processes with diverging characteristic time. We consid- 
ered two different mechanism for perturbing the system. 
For the first one, which we termed phenomenological, we 
rederived the results obtained in Refs. [16, 17] and a Fluc- 
tuation Dissipation Relation equal to the one that can be 
demonstrated to be valid for Hamiltonian systems per- 
turbed out of canonical equilibrium [14]. For the second 
approach, which we termed dynamical and consider more 
realistic because it involved a perturbation of the statis- 
tics of the "events" generating the random walk process, 
we derived an exact analytical expression for the response, 
and showed that a different FDR, Eq. (6) is fulfilled. 

Concluding Remarks. — We devote these conclud- 
ing remarks to discussing the possible consequence of the 
main result of this article: Eq. (6). The first important re- 
sult is that the equivalence between GME and CTRW not 
only is questionable, but is incorrect. This is expected to 
have significant consequences on the foundation of a the- 
ory for BQD phenomenon. Altough the response to exter- 
nal perturbation of single quantum dots within a Hamilto- 
nian model has been addressed in several works (e.g. [19]), 
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this has not been achieved in the BQD "regime", which 
occurs for quantum dots embedded in a disordered envi- 
ronment. Our result indicates that the derivation of the 
BQD phenomenon from a Hamiltonian picture is not only 
difficult but impossible, thereby locating the intermittent 
fluorescence at the level of those emergent cooperative pro- 
cesses that determine rcductionism failure. 

We speculate here also a second possible consequence of 
the main result of this article. This is connected with the 
concept of an effective temperature such as was introduced 
into the physics of glassy systems. In 1997 the authors of 
Refs. [20,21], on the basis of earlier work [22] (see also [23]) 
argued that it is convenient to replace Eq. (3) with 

RAB{t,t')^(3K{t,t')-^{A{t)B{t')) t>t', (40) 
at' 

where the function K{t,t') has the role of defining an ef- 
fective temperature depending on the age of the system 
t'. The effective temperature has then the role of link- 
ing the out-of-equilibrium configurations occupied by the 
system during its relaxation at a given temperature, to 
equilibrium configurations relative to a different tempera- 
ture. There has been an intense research activity in this 
direction, and we limit ourselves to quoting some papers 
representative of the experimental and theoretical work 
done in this direction [24-27], with controversial results 
including the report of no deviation from the FDT over 
several decade in frequency [26]. Is the generalization of 
FDT derived in this Letter with exact dynamic arguments 
compatible with the concept of an effective temperature as 
introduced in Eq. (40)? Consider Eqs. (15) and (16) from 
which the following relations with the survival probability: 



d^(MO 

dt 
d^{t,t') 

dt' 



yj{t,t') 
P(t')^{t-t'), 



(41) 



are determined. Thus, the result of this letter can be 
expressed in the form of Eq. (40) (/3=1, A = B = ^) 
with 

Kit,t')^i;it,t')[P{t')^5j{t-t')]-\ (42) 

The introduction of an effective temperature in this case 
has the role of linking our description to one compatible 
with the ordinary FDT obtained in Refs. [15,16], which in 
turn is directly linked to a Hamiltonian derivation. 
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